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Abstract 

Electronic properties of disordered binary alloys are studied via the calcula- 
tion of the average Density of States (DOS) in two and three dimensions. We 
propose a new approximate scheme that allows for the inclusion of local order 
effects in finite geometries and extrapolates the behavior of infinite systems 
following finite-size scaling ideas. We particularly investigate the limit of the 
Quantum Site Percolation regime described by a tight-binding Hamiltonian. 
This limit was chosen to probe the role of short range order (SRO) properties 
under extreme conditions. The method is numerically highly efficient and 
asymptotically exact in important limits, predicting the correct DOS struc- 
ture as a function of the SRO parameters. Magnetic field effects can also be 
included in our model to study the interplay of local order and the shifted 
quantum interference driven by the field. The average DOS is highly sensitive 
to changes in the SRO properties, and striking effects are observed when a 
magnetic field is applied near the segregated regime. The new effects observed 
are twofold: there is a reduction of the band width and the formation of a 
gap in the middle of the band, both as a consequence of destructive interfer- 
ence of electronic paths and the loss of coherence for particular values of the 
magnetic field. The above phenomena are periodic in the magnetic flux. For 
other limits that imply strong localization, the magnetic field produces minor 
changes in the structure of the average DOS. 
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I. INTRODUCTION 



Considerable efforts have been made in the understanding of localization properties of 
electronic states of disordered systems since the work by Anderson.0 Being the simplest one 
to investigate the above phenomena, the Anderson model is formulated through a Tight 
Binding Hamiltonian with diagonal disorder and hopping t to nearest neighbors. The site 
energies are randomly distributed within some range W, and W/t is the relevant parameter 
to distinguish between the strong and weak localization regimes. Anderson showed that 
disorder, in three dimensional systems, drives a metal-insulator transition, with a mobility 
edge separating localized from extended electronic states. 

The one parameter scaling theory, introduced by Abrahams et a/.,i predicts that all states 
are localized for dimensions smaller than or equal to two. Dimension two then appears as 
a marginal case. The one-parameter scaling hypothesis assumes that the conductance of a 
finite sample scales in a universal way, and has been supported by several experiments done 
in the early 1980s,! and by a large number of numerical works. However, it has been 
suggested that the theory breaks down when electron-electron interactions are included, or 
when an external magnetic field is applied.EB Remarkable new experiments on 2D electron 
systems,^ in samples with much higher electron mobilities than the previous ones, do show 
a clear evidence for a metal-insulator transition at low electron densities, in the absence of a 
magnetic field. The effect has been ascribed to the dominant role of Coulomb interactions, 
and a new scaling theory has been proposed for interacting electrons,^ leaving open the 
possibility of a 2D metal. 

Within the one-electron scheme, it was well known that in 2D, the effect of back scattering 
from impurities is enhanced, leading to an insulating behavior at zero temperature. Electron 
trajectories interfere constructively in the backward direction, and one gets localization for 
any degree of disorder. Transport properties with an external magnetic field are useful 
probes to test the above effects. Quantum interference of electron trajectories, modified by 
the applied field, may change the localization properties of the electron wave function. One 
evidence of this fact is the observation of magnetoresistance oscillations with an applied 
field in thin disordered metallic cylinders.0'0 With this geometry, and orienting the field 
along the cylinder axis, one can select the phase of the electron orbits in order to observe 
the phenomenon macroscopically. For two-dimensional systems, like the ones realized in 
semiconductor heterostructures, it has been observed a metal-insulator transition induced 
by the magnetic field.li! Samples that, at sufficiently low temperatures and zero field, exhibits 
a localized electronic behavior, present evidence of extended states at strong magnetic fields. 
Other evidences for a transition to a metallic regime have been observed in experimental 
studies of the magnetoresistance and magnetoconductance.BEl As a general rule, quantum 
interference effects are invoked to explain the derealization of the electronic wave function 
with the field.l The reentrant behavior of the magnetoresistance found by some authors 
at very high magnetic fields ,00 is associated with the shrinkage of the electronic wave 
function. Ha 

One effect marginally taken into account, and that may have important consequences, 
even within the one-electron picture, is the inclusion of Short-Range Order (SRO) correla- 
tions in the electronic properties of disordered systems. Real disordered systems are not 
necessarily random, but may present different degrees of short-range order. In a binary alloy, 
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for given concentrations of the atomic species, one may get configurations ranging from seg- 
regation to complete order as a function of the SRO properties. In this work we investigate 
the role of the above effects and its interplay with the action of an applied magnetic field, in 
The Quantum Site Percolation limit. The latter regime corresponds to the case where the 
two site energies, for a binary alloy, are very far apart in the energy scale, and percolation 
of the wave function is only possible through atoms of the same species. Thus, this limit 
favors localization when short range correlates atoms of different kind as nearest neighbors. 
We remark that this case is difficult to calculate within other approximations of nonlocal 
character .0 We calculate the average electronic Density of States (DOS) using a previously 
developed method reported in Ref. |19[ Our approach relies on one dimensional numerical 
techniques, and is efficiently applied to long stripes and bars of finite cross section with 
the idea that two and three dimensional lattices may be obtained by joining infinite linear 
chains. From grounds based on Statistical Mechanics results,^ this geometry, generically 
called as cylindric, is expected to converge more rapidly to the bulk limit than the simpler 
block geometry. Indeed, for the disordered case, the bulk regime can be extrapolated with a 
relatively small number of chains. SRO effects are included following the Kikuchi method,il 
where correlations are systematically handled. The last step in the calculation involves a 
configurational average, where the probability of a given configuration is obtained in terms 
of the short and long range order parameters. The electronic DOS calculated through the 
above procedure, in spite of being an average, is rich in structures, and as we will show 
below, is highly susceptible to changes in local properties. 

Magnetic field effects are considered with the introduction of a phase in the electronic 
hopping term! via the so called Peierls substitution.il This phase is related to the magnetic 
flux across an elementary cell, as given by the well known Aharonov-Bohm effect .Ej The 
lattice model approach adopted in this paper is known to accurately describe the lowest 
Landau level in the situation where the disorder broadening is such that no overlap of mag- 
netic sub-bands occurs. We are also neglecting electron-electron interactions and electron 
spin (spin-split Landau levels), and all the effects found in the present contribution are ex- 
clusively due to disorder in a one-particle lattice model. We found that magnetic field effects 
are stronger near the segregation limit, where one expects weak localization properties and 
a higher degree of coherence in the electronic wave function. An oscillatory behavior of the 
electronic DOS is observed as a function of the field, with surprising effects on the structure 
of the DOS. 

Our paper is organized as follows: In the next Section we present the method to calcu- 
late the DOS. As illustrative examples, we show results for one-component non-disordered 
samples in two and three dimensions. The two-dimensional disordered case including SRO 
effects is also shown for completeness. In Section III we extend our calculation to three- 
dimensional disordered systems in the presence of an applied magnetic field. Finally, in 
Section IV we present and discuss numerical results for several examples of the latter. 

II. A NEW APPROACH TO CALCULATE THE DOS 

The computational scheme that will be employed here, was developed previously by 
the authors to study the electronic DOS in two-dimensional disordered binary alloys.113 For 
completeness, we discuss in this Section some details of the method and present the extension 
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to three dimensions in the next Section. 

A tight-binding Hamiltonian for a single s-atomic-like orbital, in a hypercubic lattice 
structure in d dimensions, yields the following dispersion relation in fc-space: 



E (k) = E + 2 V cos ha , (1) 



i=i 



where E is the atomic energy, V is the hopping (kinetic) term for nearest neighbors in a 
pure system, and a is the lattice parameter. The DOS for d- dimensional systems can be 
obtained recursively from the DOS in (d — 1) dimensions, in the form of a convolution with 
the DOS of a linear chain 

V {d \E) = J dE'V^ d - l \E') V {l \E - E'). (2) 

This result is shown analytically in the Appendix for two dimensions. For a finite system 
with the geometry of a hyper- cylinder, like a strip with finite width, the DOS is obtained 
as a discrete convolution of the DOS of the linear chain with the distribution of eigenvalues 
for the finite transverse section of the cylinder 

V { Z mder {E) = £ vtl n {Ei) V^\E - Ek) (3) 

i 

In finite systems, strong oscillations are present in the DOS, specially near the points where 
van Hove singularities are developed in the infinite limit (in analogy with the so called 
Gibbs phenomenon for Fourier series). The above oscillations are typical size effects, with 
the amplitude being damped when the size of the finite section is increased.@ In Fig.|I] we 
display the DOS for the square lattice for wide but finite width stripes, along with the 
infinite two dimensional case (uppermost figure). The corresponding three dimensional case 
(cubic lattice) is depicted in Fig.|2|. We compare finite bar geometries with the infinite three 
dimensional limit. The energy scale is given in units of V, the hopping integral. The large 
amplitude of oscillations are due to coherence effects in pure samples, and as we will show 
later, size effects decay more rapidly in disordered systems. In addition, as remarked in the 
Introduction, the cylindric geometry, with an infinite dimension along the axis, is expected 
to converge considerably faster to the bulk limit than the block geometry. This latter fact 
allows the use of relatively small samples to extrapolate the DOS of disordered alloys. One 
comment is in order here. As shown in expression (^), the DOS for a finite system is obtained 
as a weighted superposition of one-dimensional densities. The interchain coupling displaces 
the center of gravity of the one-dimensional DOS to the renormalized value E^ which takes 
into account the true coordination number and the degrees of freedom of the transverse 
section of the system. The weight in fl3|) for each chain is exactly the discrete density of the 
eigenvalues {Ei}. 

The above exact analytic results are extended heuristicahV to treat disordered systems. 
In this Section we show the procedure for two dimensions,!^ where the transverse part of 
the system is a finite chain. Local order effects are handled with a single variable, the so 
called Cowley's parameter 



*=l-m, (4) 

CACB 
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with pab being the probability of having an AB pair, and ca and cb the concentrations 
for species A and B respectively. With the above definition, a ranges from —1 (limit of 
complete local order) to 1 (limit of segregation), passing through intermediate cases with 
variable local order that include the random limit for o = 0. For Markovian chains, SRO 
effects were embodied in a numerically exact algorithm based on the so called Negative 
Factor Counting (NFC) method!! to calculate the DOS for such systems. Our scheme thus 
relies on exact calculations for D^(E;a), the DOS of disordered infinite chains with local 
order described by a. To get the analogous of expression (|3|) for a disordered two-dimensional 
stripe of width M, we generate all the configurations of a finite chain of M sites, and get 
the corresponding eigenvalues and their multiplicities. For a given configuration, the DOS 
is obtained as 

M 

vfl ip (E; M,c) = J2 d^ ctlon (E i; M, c) V^(E - E t ; a), (5) 

i=l 

where the V®(E-Ei\ a) is an infinite chain DOS with the center of gravity displaced to E i: 
the eigenvalues of the finite transverse chain for the given configuration. Each term of the 
sum (||D is weighted by dg ed . ian {Ei] M,c), the discrete distribution of the eigenvalues {Ei}. 
The final DOS is an average over all the possible configurations: 

V ( 2 ip (E) = £ P c V&tpiE; M, c) , (6) 

c 

with P c equal to the probability of a given configuration c. The latter is obtained, for 
Markov chainsj^l as a weighted product of pair probabilities and depends only on a and 
the concentrations. In Fig.|3| we show the lowest sub-band for 50-50% alloying with different 
degrees of SRO in two dimensional strips of width M = 20. In spite of the finite size, the 
three cases already display two dimensional features characteristic of the infinite limit. The 
DOS structure changes dramatically as a function of the local order, with highly localized 
states near the random and ordered limits. The upper figure is close to the segregation limit 
and resembles the pure case shown in Fig.[l|. For this regime, the electronic wave function 
may percolate through atoms of the same type in very large islands. In contrast with one 
dimension, the above properties are typical of the two dimensional topology, with many 
different percolation paths for extended electron trajectories. 

For three dimensions, the finite section is of two dimensional character and SRO proper- 
ties are handled with more parameters that include the possibility of closed path correlations. 
This is done using the Kikuchi method^ and is described in the next Section. 



III. THE INCLUSION OF MAGNETIC FIELDS AND SRO PROPERTIES IN 
THREE-DIMENSIONAL GEOMETRIES 

Following ideas developed in Finite-Size Scaling methods,! our system is taken as a bar 
in three dimensions, with finite dimensions in the x and y directions, and infinite extension 
along the z direction. Periodic boundary conditions are imposed for the finite section of the 
system. The magnetic field is applied parallel to the z axis. 

We model our system using an Alloy Tight Binding (ATB) Hamiltonian: 
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n = J2 E i l*> V\ + E v ^ 3 N) 01 + ^ c -> (7) 

<i> <ij> 

where off-diagonal elements are taken different from zero only for nearest neighbors, and faj 
is a phase that represents the effect of the magnetic field in an elementary cell with the lattice 
parameter as a side (Peierls substitution^). For simplicity, we choose the hopping amplitude 
as real along the x axis (0 = 0), and complex along the y axis. The phase difference between 
two successive bonds is related to the magnetic flux enclosed in an elementary cell, 

(f>n+l - 0n = Y c J ^ ' ^ ' ( 8 ) 

where A is the associated electromagnetic vector potential. Periodic boundary conditions 
for the finite section of the system induce a quantization of the magnetic flux that depends 
on the size of the sample. One can get rid of the above constraint by choosing free ends, 
but in this case finite size effects decay more slowly. 

As mentioned in the previous Section, we studied the limit corresponding to the Quantum 
Site Percolation case, where \Ea — Eb\ » V. SRO is treated using the Kikuchi method 
with closed square diagrams as highest order correlations.il This approximation is known 
to correctly describe the two dimensional topology. For each square we have 16 different 
configurations. These configurations are displayed in Fig.[|, with the respective probabilities 
(^s) and degeneracies (/3 s). 

Normalization conditions require: 

C A + C B = 1 (9) 

for concentrations, 



'■i 



Paa + Pab = c A 

Pbb+Pab = c b (10) 



for pairs, and 



Paa = zi + 2z 2 + z 3 

Pab = z 2 + z 3 + + z 5 (11) 
Pbb = z 3 + 2z 5 + z 6 

for squares. 

In all, we get 5 independent parameters, which are chosen to be Pab, z>3, %4, £i and £ 2) 
with the last two defined as: 

Ci = c a ~ c B 

6 = *2 - z 5 . (12) 

While the first three parameters (pab, z%, z±) are related to correlations of short range, the 
last two (|12D handle long range order. 
Dependent variables are written as: 
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ca = 


(l + eO/2 


c B = 


(l-Ci)/2 


Paa = 


(l + ^-2p AB )/2 


Pbb = 


(i-a-2 PAB )/2 


Zl = 


(1 + & - 4 Pi4 B + 2z 4 - 2£ 2 )/2 


z-i = 


(pAB-^3-^4 + e2)/2 


z 5 = 


(pas - 23 - z± - 6)/ 2 


ZQ = 


(1 - & - 4^ + 2^ 4 + 2£ 2 )/2 



Within the approximation adopted, a given configuration of the finite section is obtained 
by joining squares. The corresponding probability is obtained as a product of conditional 
square probabilities, where the constraint is that two squares are linked by a common bond. 
This approach has been devised in analogy to Markovian chains,@ but with the square as 
the basic correlation. 

As it will be seen in the next Section, our method allows to include local properties under 
the coupled effects of variable disorder and application of a magnetic field. 



IV. NUMERICAL RESULTS AND DISCUSSION 

In the examples that follow, we consider 50-50% alloying and several degrees of local 
order. The site energies are taken to be Ea = and Eb = 1000, in units of the hopping 
term V. In all the figures, we only display the lowest sub-band, the one centered at Ea- 
Fig.|5] shows examples calculated in the absence of a magnetic field: (a) near the ordered 
limit; (b) random disorder; and (c) near the limit of segregation. The behavior is similar to 
the one shown previously for two dimensions. Near the segregated limit the DOS for each 
sub-band resembles the one of the pure case, suggesting the presence of extended electronic 
states. On the other hand, when the system is near the ordered limit, the DOS shows the 
typical structure of localized states, as it should be expected for two kinds of atoms with 
very different site energies. In this instance, the probability for one state to propagate is 
nearly zero, due to almost infinite barriers. The authors are currently calculating other 
cases where the centroids approach to each other, and the bands overlap. In this situation, 
quantum diffusion is possible by means of tunneling through high but finite barriers, and 
a derealization onset is expected as a function of the band separation.il To illustrate the 
role of local order, in Fig.|6| we display the DOS values for three nearby energy levels as a 
function of the correlation a, in two dimensions. One of the levels corresponds to a sharp 
peak in the middle of the band in the lowest part of Fig.[| for a = 0. 

Size effects are pictured in Fig.[7] for the random alloy. It is worth noting that the DOS 
of finite systems already present the whole structure of the infinite case, even for small 
transverse sections. This seems to be characteristic of the bar geometry, where one of the 
dimensions is infinite and the bulk limit is extrapolated with not too big cross sections. 
Due to numerical limitations, the next examples in the presence of a magnetic field, are 
calculated in finite geometries only (sections of 12 atoms). 

We now analyze the role of the magnetic field. In Fig.|| we show the DOS for a case 
near the segregation limit. An oscillatory behavior of the DOS with the field is observed, 
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and the figure depicts cases up to half the period in the flux. The values are adjusted to 
periodic boundary conditions. The effect is paramount at the center of the sub-band, with 
the appearance of a gap for A0 = tt. The phenomenon is due to destructive interference 
of the electronic paths in the presence of the field. The band width is also affected, with a 
striking reduction at half a period. 

On the other hand, tiny effects are observed near the random and ordered limits, with 
no sensitive changes in the structure of the DOS as a function of the field. Other cases with 
different concentrations were also calculated. As a general rule, the DOS structure is more 
affected by the field near the so called weak localization regime, when the concentration 
of one of the species is small. Localization and derealization of electronic states may be 
interpreted as pure quantum effects due to destructive or constructive interference of the 
different electronic paths. This interference is shifted by the application of a field. In Fig.|9|, 
we show the effect of the magnetic field on a particular level as a function of short-range 
correlations. The figure displays a trajectory in our parameter space, ranging from the local 
ordered limit to the segregated case. This trajectory, for short, is labeled with the symbol 
a. The figure pictures in a different way our comments above: the DOS oscillates as a 
function of the field, but the effect predominates in the vicinity of the segregated limit (for 
0.5 < a < 1 in the figure). 

In summary, it has been shown that the method proposed is numerically efficient, repro- 
duces the DOS in important limits, and predicts new behaviors as function of the SRO and 
magnetic field. To theoretically investigate localization properties, one needs to calculate 
the electron wave function, or get indirect information via the so called localization length. 
In two and three dimensions, most of the works in the literature concentrate around simula- 
tions or other numerical techniques. As a common characteristic, they are based on Finite 
Size Scaling methods, but most of the cases investigated correspond to random disordered 
systems, represented generally by the Anderson Model or the Quantum Percolation Modell 
Our results, in spite that we calculate the average DOS, suggest that SRO effects may play 
an important role in the localization of electronic states. A phase diagram in which the 
degree of localization is displayed, may depend not only on the concentrations, but also on 
the SRO parameters. In order to illustrate this issue, we suggest numerical computations of 
the localization length in systems like the ones considered here, with equal concentrations 
of the species, very different site energies, and variable local order properties. Changes in 
the localization properties will then be entirely ascribed to short-range correlations. 

The applied magnetic field strongly changes the electronic properties in the weak localiza- 
tion limit. On the other hand, in the opposite regime, the effects of the field on the average 
DOS are negligible. Huge magnetic fields are involved in the above phenomena due to the 
small size of the elementary cell. However, the effects predicted here may be brought to a 
range susceptible of observation in the case of superlattice structures, or in geometries where 
the electronic paths enclose macroscopic areas that contribute to the flux. We note that the 
present state of the art in nanostructure technology allows the fabrication of modulated 2D 
systems with lattice parameters of the order of 100 nm or more (see for instance Ref. 
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APPENDIX A: EXACT CALCULATION FOR THE DOS IN TWO DIMENSIONS 

The DOS for the pure cubic lattice in d dimensions was given in Section II by eq.(^). 
We show explicitly this result for the two-dimensional case. 

If we consider that states {|m, n >} represent the basis in the site representation, where 
the indexes are integers along the x and y axes respectively, the Hamiltonian reads: 



H = E mn \m, n) (m, n\ + V] t x \m, n) (m, n + 1| + h.c. 

<n> 

+ t y \m, n) (m + 1, n| + h.c. (Al) 



<in,n> <n> 



<m> 



Taking the Fourier transform in relation to the x direction: 



\m,n) = -= e ikn \m,k) , (A2) 



N <k0 > 



we get: 



TC = )XE Q + 2t x cos k) \m, k) (m, k\ + 

m,k 

Y t y (\m, k) (m + l,k\ + \m + 1, k) (m, k\) . (A3) 



m,k 

The Green's function in the basis {\m, k >} is given by:0 

= 1^ k ) (A4) 

mk z-E{k, m y 

where E(k, m) are the eigenvalues: 

E(k,m) = E + 2t x cos k + X m = E x (k) + A m , (A5) 

with: 

E x (k) = E + 2t x cos k, (A6) 

and A m corresponding to the eigenvalue of an effective chain with site energy E + 2t x cos k 
and hopping t y between first neighbors along the y direction. 
The Green's function can be rewritten as: 
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G{z) = J2 |m ' fc)(m,fc| (A7) 

m,k Z ~ E x {k) — X m 

with the diagonal element: 

L r 1 

(m, n\ G(z) |m, n) — — — / dk — —r (A8) 

x ' 1 v ; 1 ' ' 2irN Jbz z-X m - E x (k) v ' 

The DOS and the Green's function satisfy the relational: 

V(E) = —ImY (m,n\ G(E) \m,n) , (A9) 

7T — ' 

" m,n 

with: 

-t^ttY (m,n\G(z)\m,n) = ^~Y — T dk — — , (A10) 

NM ^ n v ' 1 y 1 1 ' 1 M^27iJ-n z-X m -E -2t x cos{k) ' v ; 

and subsequently the DOS is given by: 

V ( E ) = ^E n ( A "0 V ldim (E-\ m ). (All) 

where T>i<ii m {E — X m ) is the DOS of an infinite linear chain centered at X m , and n(X m ) is 
the degeneracy of the eigenvalue. 
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FIGURES 



FIG. 1. DOS, in arbitrary units, for a pure two-dimensional system: (a) upper figure, the exact 
two dimensional infinite sample; (b) middle, a strip with the finite transverse length M = 500; (c) 
lower, a strip with the finite transverse length M = 100. In all the cases, we take Ea = Eb = 
and Vij = 1, with energy steps equal to 5.10~ 4 . 

FIG. 2. DOS, in arbitrary units, for the pure three-dimensional system: (a) upper figure, 
the exact three dimensional case obtained as the convolution of a two-dimensional DOS with a 
one-dimensional one; (b) middle, bar with N x M = 48; (c) lower, a bar with N x M = 24. For 
all the examples, Ea = Eb = and = 1, with energy steps equal to 5.10~ 4 . 

FIG. 3. DOS, in arbitrary units, for the square lattice with several degrees of SRO. All the cases 
shown in the figure were calculated for infinite stripes of width M = 20 chains. Concentrations 
are equal for both species and the Cowley parameter is given in each figure. Other parameters 
are Ea = 0, Eg = 1000, V = 1, with energy steps equal to 5.10~ 2 . We only display the lower 
sub-band: (a) upper figure, a = 0.8, tendency to segregation; (b) middle, a = 0, random alloy; (c) 
lower, a = —0.8, tendency to ordering. 

FIG. 4. Possible square configurations for two kinds of atoms (filled and empty circles). We 
denote by Z{'s the probabilities of each configuration, and the numbers indicated in the last column 
are the respective degeneracies. 

FIG. 5. DOS, in arbitrary units, for the three-dimensional disordered case with several degrees 
of SRO. All the examples correspond to the infinite system, and were obtained by the convolution 
technique explained in the text. The concentrations of both species are equal, and as in previous 
examples, we use Ea = 0, Eb = 1000, V = 1, with energy steps equal to 5.10~ 2 . We only 
display the lower sub-band for the cases: (a) tendency to order; (b) random alloy; (c) tendency to 
segregation. 

FIG. 6. DOS, in log scale, for the two dimensional case, for three nearby energy points, as a 
function of the short-range parameter a. 

FIG. 7. Size effects for the random disordered binary alloy for the three dimensional case 
(case (b) of Fig.||). For finite systems we use the following correlation parameters: pab = 1/4, 
Z3 = Z4 = 1/16, and £1 = £2 = 0. Note that the lower DOS, for a bar of cross section of 12 atoms, 
already presents, except for the amplitude of some peaks, the main structure of the infinite system 
density. 
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FIG. 8. DOS, in arbitrary units, for a three-dimensional bar (N x M = 12) for different values 
of the applied magnetic field. The case depicted corresponds to a disordered binary alloy near the 
segregation limit for equal concentrations of both species. Correlation parameters are pab = 0.05, 
zs = Z4 = 0, and £i = £2 = 0. We only show the lower sub-band and the magnetic flux is indicated 
in each figure. Due to periodic boundary conditions for the finite section of the bar, the magnetic 
phase factor is quantized in units of ir/3. 

FIG. 9. DOS, in arbitrary units, at a particular energy point (E = 3.0), for a three-dimensional 
bar (N x M = 12), for different values of the applied magnetic field, as a function of the short-range 
order. 
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